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ABSTRACT 



Aims. We investigate the emission mechanism and evolution of pulsars that are associated with supernova remnants. 

Methods. We used imaging techniques in both the optical and near infrared, using images with very good seeing (<076) to study 

the immediate surroundings of the Crab pulsar In the case of the infrared, we took two data sets with a time window of 75 days, to 

check for variability in the inner part of the Crab nebula. We also measure the spectral indices of all these wisps, the nearby knot, and 

the interwisp medium, using our optical and infrared data. We then compared the observational results with the existing theoretical 

models. 

Results. We report variability in the three nearby wisps located to the northwest of the pulsar and also in a nearby anvil wisp in terms 
of their structure, position, and emissivity within the time window of 75 days. All the wisps display red spectra with similar spectral 
indices (a v = -0.58 + 0.08, a v = -0.63 + 0.07, a v = -0.53 + 0.08) for the northwest triplet. The anvil wisp (anvil wisp 1) has a 
spectral index of a v = -0.62 + 0.10. Similarly, the interwisp medium regions also show red spectra similar to those of the wisps, with 
the spectral index being a v = -0.61 + 0.08, a v = -0.50 + 0.10, while the third interwisp region has a flatter spectrum with spectral 
a v = -0.49 + 0.10. The inner knot has a spectral index of a v - -0.63 + 0.02. Also, based on archival HST data and our IR data, we 
find that the inner knot remains stationary for a time period of 13.5 years. The projected average velocity relative to the pulsar for this 
period is < 8 km s _1 . 

Conclusions. By comparing the spectral indices of the structures in the inner Crab with the current theoretical models, we find that 
the Del Zanna et al. (2006) model for the synchrotron emission fits our observations, although the spectral index is at the flatter end 
of their modelled spectra. 

Key words. Supernova remnant - Crab - pulsar - Optical - Infrared 



1. Introduction 

The Crab nebula is one of the most studied targets in the sky. It 
is a remnant from a supernova explosion that occurred in 1054. 
The Crab pulsar (my- 16), the first pulsar with optically detected 
pulses (Cocke et al. 1969), is located at the centre of the nebula, 
and is responsible for powering the nebula. The pulsar, together 
with the filaments and the continuum emitting pulsar- wind neb- 
ula (PWN), form the Crab nebula. The Crab nebula was the first 
detected astrophysical source of synchrotron radiation, as sug- 
gested by Shklovsky (1953) and was confirmed by polarisation 
observations a year later by Dombrovsky (1954). The inner part 
of the Crab nebula is a region consisting of jets, a torus of X-ray 
emission (Aschenbach & Brinkmann 1975), and complexes of 
sharp wisps (Scargle 1969). The structures give a unique oppor- 
tunity to study the pulsar, along with its wind and the interaction 
it has with the rest of the remnant. 
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Copenhagen. 



Evidence of activity in the Crab pulsar and its surroundings 
has existed for a few decades. Variability in the structures around 
the pulsar were mentioned by Lampland (1921) and by Oort & 
Walraven (1956). A very comprehensive study has been done by 
Scargle (1969), who found that the wisps may show relativistic 
motion, but conversely always seem to be more or less located 
at the same position. Using the WFPC2 onboard HST high res- 
olution surveys were conducted by Hester et al. (1995; 2002). 
These data were complemented with X-ray data from ROSAT 
and Chandra to get a complete picture of the Crab pulsar and its 
immediate surroundings. The most outstanding discovery is the 
presence of two knots that are situated 0765 and 378 southeast of 
the pulsar. The X-ray data have also shown many knots around 
the pulsar that show variability with time. For a complete look at 
the Crab pulsar and its nebula see Hester (2008). 

Sollerman (2003) used ISAAC on the VLT, to conduct an 
IR survey on JHK S bands in order to investigate the Crab pulsar 
wind nebula further. Using the VLT data, together with archive 
HST optical data he measured a spectral index for the knot equal 
to a v = -0.8. Melatos et al. (2005) went further into investigat- 
ing the Crab PWN using IR images to check for short time- scale 
variations around the Crab pulsar. They found variations in the 
emission of the Crab wisps on a time scale of 1.2 ksec of ±24% 
± 4% in the K band and ±14 ±5% in the J band. Apart from the 
variation estimations, Melatos et al. (2005) calculated the spec- 
tral indices from the structures near the pulsar, finding values for 
the spectral index that vary from a v = -0.16 for the rod (this the 
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brightest feature in the shock located ~ 4" in a southeastern ori- 
entation from the pulsar) down to a v = -0.76 for the innermost 
knot. 

Despite the amount of work that has been done on the Crab, 
there are still many unanswered questions regarding its nature 
and, in general, about the supernova explosion mechanism and 
how the pulsar emission mechanism operates and evolves. Also 
there is not any model that can fully explain the observed fea- 
tures of the emitted radiation, especially in the inner part of the 
PWN. In this article we present our near infrared and optical 
survey that we conducted using the Nordic Optical Telescope in 
La Palma, Spain, with a gap of two and a half months (in the 
case of the NIR) in order to check for variability in the enviro- 
ment around the Crab pulsar. In Sect. 2 we present our observing 
procedure followed by the data reduction process. In Sect. 3 we 
present our photometric results, and in Sect. 4 we present argu- 
ments for the nature of the wisps and the knots. We make our 
conclusions in Sect. 5. 




Fig. 1. 1-band image of the Crab nebula taken with ALFOSC at 
the NOT on 7 December 2007. The pulsar has been removed 
with the PSF to reveal the nearby knot that lies at a distance 
of Of! 65. The field of view is ~ 45x45 arcseconds, with north 
pointing upwards and east to the left. The areas that we used to 
measure the emission from the wisps and the interwisp regions 
are marked with circles. 



2. Observations 

The NIR observations were carried out on two seperate occa- 
sions. The first took place on 30 September 2007 and the sec- 
ond on 13 December 2007. Both runs were done using the high- 
resolution camera of NOTCarrQon the Nordic Optical Telescope 
(NOT) in Observatorio del Roche de los Muchachos in La Palma 
Spain. The observations used an H-band filter (A c = 1.626 yum, 
FWHM = 0.296 yum) and a filter in the K s band (A c = 2.140 
fjm FWHM = 0.310 /mi). The detector of the NOTCam is an 
18.5 fimx 1024 x 1024 HgCdTe (Hawaii) array. The pixel scale 

1 http://www.not.iac.es/instruments/notcam/ 




Fig. 2. Upper panel: H-band images of the inner part of the 
Crab nebula taken with the high-resolution camera of NOTCam 
(07078/pix) at the NOT on 30 September and 13 December 
2007. The pulsar appears to be elongated due to the presence of 
the nearby (0765) knot. The wisps that surround the pulsar are 
prominent on both frames. By setting the cuts in an appropriate 
manner, the knot can also be visible in these images. The knot 
can be clearly seen after the subtraction of the pulsar on the im- 
age from 30 September (middle). On the right side of the panel, 
we show the H band image from the night 13 December, where 
the changes in the wisp triplet are obvious in terms of brightness 
and structure especially for the outer wisp. Similar changes can 
be seen in the anvil region too. (See Fig. 1 for a description of the 
Crab nomenclature) Lower panel: K s band images of the inner 
part of the Crab nebula from the same nights. On the left side of 
the panel, we show the image from the night of 30 September. In 
the middle of the panel we show in the inner knot after we sub- 
tracted the pulsar with the PSF, and finally in the right side the 
image from the night of 13 December. In all images, the field of 
view is ~ 35x32 arcseconds, and north is up, east is to the left. 



of the HR camera is 07078/pix which gives a field of view of 
80" x 80" . The first observing run used the Engineering Grade 
Array, while the second run used the new Science Grade Array 
for NOTCam. 

The observations were performed in the same manner for 
both NOTCam runs. Because the whole field of view of the HR 
camera contains extended emission from the nebula, we used 
beam- switching between target and sky, located 180" towards 
the south, in order to account for a proper background sub- 
traction. Each position (ON target, as well as OFF target) was 
exposed for 80 seconds using a ramp-sampling readout mode 
where the array is read nondestructively every 10 seconds. A 
linear regression analysis is performed on the 8 readouts to give 
the final image. 

We achieved a total exposure time of 1500 sec for the H band 
and 750 sec in the K s band, partly limited by weather on the first 
night. The average seeing on the night of 30 September was 0745 
and the data were taken under photometric conditions. On the 
second night, due to variable seeing (0750 < seeing < 0780), we 
selected the sharpest images only (< 0765), giving an effective 
total exposure time of 640 sec in the H band and 800 sec in the 
K s band. The seeing in the coadded IR images on the night of 30 
September is ~ 0745, while for the night of 13 December is ~ 
0755. On both nights we observed the photometric standard star 
AS05 (Huntet al. 1998). 

Optical observations were carried out on 7 December 2007 
(I band) and 20 December 2007 (U & V) using the Andalucia 
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Source Filter Exp. time Eff. Exp. Time Seeing Date 
(Seconds) (Seconds) (Arcsec) 



Source 


Date 


H-band 


K s -band 


Ml 


2007/09/30 


14.14 + 0.02 


13.64 + 0.05 


Ml 


2007/12/13 


14.10 + 0.02 


13.58 + 0.05 



The photometric zeropoints of the two infrared filters were 
obtained from the photometric standard AS05 (Hunt et al. 1998), 
and they are found to be 23.70 ±0.01 mag. for the H band and 
23.16 ±0.01 mag. for the band for the data that we got on 30 
of September and 24.10 ±0.01 mag. for the H band and 23.58 
±0.01 mag for the for the data from 13 December. (The dif- 
ference in the photometric zero points is due to the fact that we 
used two different arrays. Also the estimated zeropoints are for 1 
electron per second) The magnitudes of the pulsar on both occa- 
sions were calculated with the standard PSF subtraction method 
within IRAF's DAOPHOT (Stetson 1987). Our H and K s mag- 
nitudes of the pulsar give slightly lower fluxes than the 2MASS 
catalogue by 0.09 and 0.14 magnitudes in the H and band, 
respectively (see Table 2) for 30 of September and 0.05 and 0.08 
magnitudes for 13 of December. We believe that our PSF fit- 
ting has good accuracy at both epochs since the knot east of the 
pulsar is clearly visible at both epochs. The errors that we mea- 
sured for the pulsar, as given by IRAF's ALLSTAR, are lower 
than 0.05 magnitudes for both filters. This also includes the 0.02 
magnitudes uncertainty estimated from the aperture correction. 
The ALLSTAR package is the most suitable package for varying 
enviroments like the one of the Crab pulsar so no other methods 
for the PSF subtraction were tested. 

For the reductions of optical data we used THELlQ which 
is an automatic pipeline for the reduction of optical and infrared 
data (Erben et al. 2005). The data were bias- subtracted and flat- 
field corrected. Because of the dithering that we applied, we also 
made weight maps for each frame in order to monitor each pixel 
separately. We then applied an astrometric solution to our frames 
and removed the background on each frame by using the images 
themselves by subtracting a constant value. We then shifted and 
coadded the images. 

Once again, the complex enviroment around the Crab pulsar 
required constructing PSF models to measure the magnitudes of 
the pulsar in the U, V, and I bands and to isolate the knot in the I 
band. For the flux calibration of the data, we used a spectrum of 
the pulsar from Sollerman et al. (2000). For the I-band images 
from 7 December, we found a photometric zeropoint of 25.35 
mag from which we estimated the magnitude of the pulsar to be 
16.15 mag. For the U and V bands, we estimated zeropoints of 
24.90 mag and 23.75 mag, respectively, and we found the mag- 
nitudes of the pulsar to be 17.55 mag and 16.80 mag, respec- 
tively. In all three filters (U, V, and I) the estimated uncertainty 
is around 0.02 mag. We also applied the needed corrections for 
the effect of atmosperic extinction. 

3. Analysis of photometry 

3.1. The inner nebula 

Our images show most of the known features that are present 
in the inner part of the Crab nebula. In all our frames we see 
the wisps in great detail. In the case of the H band image, there 
is some contamination by the [Fe II] 1.64 fim line detected by 
Graham et al. (1990). The K s band image is dominated by the 

2 http://www.astro.uni-bonn.de/ mischa/theli.html 
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Faint Object Spectrograph and Camera (ALFOSC) at NOT0. 
ALFOSC has a pixel scale of Off 19 per pixel and it offers a field 
of view of 6 / .5x6 / .5. A series of frames were taken using an I (A c 
= 7970 A FWHM = 1570 A) continuum filter in the night of 7 
December. In total we exposed for 2900 sec in the I band in a 
13x200sec procedure plus a 300 sec exposure. The U band im- 
ages on 20 December 2007 were taken with a broad U filter (A c 
= 3550 A, FWHM = 550 A) with a total exposure time of 3600 
sec. For the images in the V band we used a Stromgren filter (A c 
= 5470 A, FWHM = 210 A) with a total exposure time of 2700 
sec. We chose these filters since, according to their profiles, there 
is no contamination from any strong emission line (e.g., Ha or 
O [III]). To reveal the nearby (i.e., ~ 076) knot we used only the 
exposures with seeing below 076. Due to seeing limitation the 
knot can not be revealed in the U and V bands images, so we 
used them just to study the wisps. For the purposes of flux cal- 
ibration of our data, we used a spectrum of the pulsar that we 
took from Sollerman et al. (2000). 

Images in the V band (A c = 5480 A, FWHM = 485 A) taken 
with the WFPC2 with total exposure time of 2000 sec, obtained 
on 14 August 1995 were requested from HST science archives 
to check that the knot next to the Crab pulsar has moved during 
the 13 -year period. 

2.1. Data reduction 

The near-IR images were reduced using IRAlQ and a set of our 
own scripts. All the images were flat-fielded by using differen- 
tial sky flats. These can be created by taking two sets of flat fields, 
one set with high counts (10,000-15000 ADU) and another with 
low counts (700-1500 ADU) if these have been obtained within 
a limited time span, which is the case for the right window dur- 
ing twilight. These frames are used to create master flats, and 
the differential approach ensures that the thermal contribution is 
subtracted. Another benefit of obtaining these series of frames 
is that they can be used to determine the bad pixels in the array. 
For cosmetic reasons we applied a bad pixel mask to remove bad 
and warm pixels from our science frames. Sky subtraction was 
done by subtracting the master sky frame that we made from the 
individual sky exposures we took. The master sky was properly 
scaled before subtracting it from each target image. Then all the 
dithered target images were flat-field corrected, aligned, shifted, 
and median combined. 



2 http://www.not.iac.es/instruments/alfosc 
1 http : //iraf . noao . edu/doc s/photom. html 
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Table 3. Proper motion of the Crab pulsar. 




Date 


Knot 


Star 1 


Star 2 


94/08/14 


07625(0.005^ 


5700 


4794 


07/09/30 


(7/635(0.012) 


5722 


5713 



Fig. 3. Top: I band image of the Crab pulsar taken with 
NOT/ALFOSC on the night of 7 December. Bottom: The same 
image after the subtraction of the pulsar with the PSF. The field 
of view is approximately 45x45 arcseconds. The residual emis- 
sion southeast of the pulsar's position is the nearby knot first 
reported by Hester et al. (1995). On both images, north is up and 
east to the left. 



synchrotron emission, but a small amount of the emission comes 
from the molecular hydrogen 1— »0 S(l) 2.122 /zm and the Bry 
2. 163 fim lines. (Use Fig. 1 to get an orientation of the photomet- 
ric analysis that we made in the inner part of the Crab nebula). 

3.2. The pulsar and the knot. 

The pulsar appears to be elongated on all images with good see- 
ing (< 076). This due to the presence of the nearby knot that 
was first identified by Hester et al. (1995) on the HST/WFPC2 
images. Hester et al. (1995) mentions the presence of a second 
knot near the pulsar, in the same orientation as the first, but lying 



closer to the anvil at a distance of 378 from the pulsar. From our 
data set there is no evidence of this knot since it is dimmer than 
the first one by a factor of almost five (Hester et al. 1995). 

After subtracting the pulsar with the PSF we fitted an aper- 
ture in order to measure the flux of the nearby knot. In the case of 
the K s band the flux of the knot is ~ 10% of the flux of the pulsar. 
This shows that in case of bad seeing conditons one would mea- 
sure a brighter pulsar. For the H band the flux of the knot is ~ 8% 
of the flux that we measured for the pulsar. The measured fluxes 
are 2.80 ± 0.15 xlO" 27 ergs sec _1 cm" 2 Hz _1 for the K s filter and 
2.65 ±0.15 xlO" 27 ergs sec^cm^Hz -1 for the H filter. Two and 
a half months later we found that the emitted flux from the knot 
is constant. The measurements from the night of 13 December 
gave fluxes of 2.90 ± 0.10 xlO" 27 ergs sec _1 cm" 2 Hz _1 for the 
K s filter and 2.75 ± 0.10 xlO" 27 ergs sec^cm^Hz -1 . Below we 
discuss our observations and results. 

The knot is not circular in shape, but has an arc-like shape. 
The latter can be used as an argument that it is not a background 
or foreground star and that it is indeed a structure associated with 
the pulsar. Another feature that supports this scenario is that this 
knot, together with the second knot that Hester et al. (1995) men- 
tioned, is aligned with the X-ray jet emanating from the pulsar. 
In the I band images, this is more pronounced since we see the 
knot aligned with an outflow of gas (optical jet) that goes through 
the anvil (see Fig. 2). Also its morphology indicates that it has 
to be a shock wave or some instability feature in the flow that is 
responsible for the jet. 

In the I band images we only managed to get two satisfac- 
tory frames in terms of seeing during the night of 7 December to 
extract the knot near the pulsar (see Fig. 3). Once again we man- 
aged to reveal the knot after the removal of the pulsar with the 
PSF. The structure of the knot appears to be the same as in the 
H and K s images. The measured distance is again found to be ~ 
0" 65, and the flux that we measured from the knot in the I band 
is 1.60 ± 0.10 xlO" 27 ergs sec _1 cm" 2 Hz _1 , while the pulsar has 
a flux of 2.45xl0" 26 ergs sec _1 cm~ 2 Hz _1 . The flux of the knot 
therefore corresponds to ~ 6.5% of the flux of the pulsar. 

We used HST/WFPC2 archival images to study a possible 
motion of the red knot with regard to the background/foreground 
stars and the pulsar. Data that date back 13.5 years from our ob- 
serving run were used for this purpose. According to the current 
measurements the Crab pulsar has proper motion of 15 mas yr" 1 
(Kaplan et al. 2008). This means that the pulsar in this amount 
of time has moved ~ Of! 2. We measured the distances between 
the pulsar and the knot at both epochs, and the two results agree 
(cf. Table 3). If the knot was a background/foreground source, 
then its proper motion would have been different than the one 
of the pulsar. This is also supported by the comparison that we 
made with stars close to the pulsar. In Table 3, we show the rel- 
ative distance between the pulsar and the knot, as well as two 
nearby stars. The pulsar-knot distance has not changed, whereas 
the Crab has moved about 072 relative to the stars, as expected. 
The distance between the knot and the pulsar that we measure 
is very close to the values estimated ~ 0765) by Hester et al. 
(1995) and Sollerman (2003). 

From our measurements we can estimate an upper limit for 
any possible relative motion in the pulsar-knot system. By us- 
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ing the relation V pm = A.lAjiD where V is transverse velocity 
in units of kms" 1 (Lyne & Lorimer 1995), one can estimate the 
shifting of the pulsar, where ji is the proper motion per year in 
mas and D is the distance in kpc. Using a distance of 2 kpc 
(Kargaltsev & Pavlov 2008) for the Crab and a proper motion 
of 15 mas/yr (Kaplan et al. 2008), the projected velocity of the 
pulsar is ~ 142 km s" 1 . The motion appears to be aligned with 
the jet axis, the pulsar moving in a northwest direction towards 
the wisp triplet. By using the uncertainties from the separation 
estimation between the pulsar and the knot, we estimate that the 
upper limit to their relative projected velocity is not more than 
~ 8 km s" 1 , which therefore makes it safe to state that the knot 
is quasistationary with regard to the pulsar, at least over a pe- 
riod of 13.5 years. Thus, assuming that the knot lies along the 
jet axis and that the jet axis is tilted -30 degrees to the plane 
of the sky (Hester 2008), as well as taking into account that the 
pulsar also moves along the jet axis, the knot trails the pulsar by 
< 10 km s" 1 , corresponding to < 6% of the space motion of 
the pulsar. 

3.3. Time variable changes in the inner nebula 

By looking at our first data set we see that there is a change 
in the region where the anvil is present (See Fig. 2 for a com- 
plete view). An obvious difference from the HST/WFPC2 V 
band image from Hester et al. (1995) is that we see two clear 
wisps in the anvil region, and we also see a presumed outflow 
of gas that is aligned with the X-ray jet (see also Fig. 3). The 
northwest direction of the pulsar shows the known features. We 
see the three wisps that are separated from each other by just 
a few arcseconds. The nearby wisp is dim on the images taken 
on 30 September. Wisp 2 is the brightest one followed by the 
outer wisp that shows some structure, where it is obvious that 
its edges are breaking. Two and a half months later, it seems 
that the enviroment around the pulsar had changed significantly. 
We cleary see that the nearby wisp west of the pulsar appears to 
have brightened and also shifted outwards. Changes are also seen 
in the wisps that are further away in the western direction. The 
brightest wisp (wisp 2) appeared to have a small change in its 
thickness, and its brightness had increased by a factor of -1.20 
in the H band, while the one that lies the furthest away, seems to 
have a single structure now. On average the flux in the H band 
is higher by a factor of ~ 1.20 on the night of 13 December. 
The flux on the K s band increased also by the same factor. The 
two wisps in the anvil region are more clearly seen now on both 
bands, and in the case of the H band image we see that the pulsar 
is surrounded by a halo, which appears to be almost 30" across. 
Other evident features in the images of 13 December are the 
outflow of gas that goes through the anvil region and a low sur- 
face brightness gas region just next to the pulsar (orientated to 
the west). In Table 4 we show all the photometric measurements 
from the wisps on both epochs and also the distance that the 
wisps have from the pulsar (See Fig. 4 to see the shift of the 
wisps between the two epochs). For the galactic extinction, we 
assumed Ry = 3.1 and E(B - V) = 0.52 using Cardelli et al. 
(1989). 

3.4. Spectral indices from the inner enviroment of the Crab 

From the ALFOSC and NOTCam (13 December) data we mea- 
sured the spectral index of the knot, the anvil (wisp 1), and the 
northwest wisp triplet, and we compared our results with the 
spectral index of the Crab pulsar and with the regions between 




Fig. 4. Stacked H band image from our two observing runs on 
30 September and 13 December. This shift corresponds to l'/2 
or 2800AU. Such a movement corresponds to a speed that is 
~ 0.2c. North is to the right and west is up. 



the wisps (interwisp medium). It is well known that the Crab pul- 
sar has a flat spectrum. The pulsar emission from y-rays to IR is 
produced by synchrotron radiation (Lyne & Graham 1998, and 
references within). Electrons following a power law N(E)dE = 
CE" r dE, where E is the energy, C a constant, and y the electron 
spectral index, will produce synchrotron radiation that obeys a 
power law F v = K(v/v ) av ergs s" 1 cm" 2 Hz" 1 . For synchrotron 
radiation, the photon spectral index a v is related to the electron 
spectral index yvma v = -(y - l)/2. 

For our measurements, we chose the same regions in all the 
images for all three wisps and we measured the surface bright- 
ness there. These regions are located on the jet axis. (See Fig. 
1 for orientation, and Tables 5 and 6 for results.) A spectrum of 
the pulsar taken from Sollerman et al. (2000) and our photomet- 
ric data were dereddened using R v = 3.1 and E(B - V) = 0.52. 
By performing a linear fit to both data sets, we found that the 
spectral index of the pulsar is a v = 0.12, while the wisps 
have spectral indices equal to a v = -0.58 ± 0.08 for wisp 1, 
a v = -0.63 ± 0.07 for wisp 2 (see Fig. 5), and wisp 3 has a 
spectral index of a v = -0.53 ± 0.08. Thus we see that the wisps 
have a rather red spectrum that is similar to the spectrum of the 
knot for which we found a spectral index of a v = -0.63 ± 0.02. 
Despite using only three points for determining the spectral in- 
dex of the knot, this resemblance can be used to give additional 
support to the pulsar-knot association, since the three wisps and 
the knot show a spectrum with similar behaviour. Finally, the 
spectral index from the anvil region (anvil wisp 1) also shows a 
red spectrum similar to those of the other wisps. We measured a 
spectral of a v = -0.62 ± 0.10. 

We also measured the flux in the regions between the wisps 
in our frames in order to compare the values with the measured 
wisp fluxes and find that these regions have a spectrum similar to 
that of the wisps and the knot (a v = -0.61 ± 0.08, a v = -0.50 ± 
0.08, a v - -0.49 ± 0.07), although the surface brightness is of 
course lower than for the wisps. 

Veron-Cetty & Woltjer (1993), who made a spectral index 
map of the whole Crab nebula in the optical by dividing it into 
small 10" x 10" boxes (see their Fig. 1), found a variation in the 
spectral index within the PWN. Their results showed an increase 
of a v as the distance increases from the pulsar, starting from val- 
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Table 4. Infrared results from the Crab PWN. 



Region 


Filter 


Flux(September 30) 


Flux(December 13) 


Distance(September 30)^ 


Distance(December 13) 


Knot 


H 


2.65(0.15]^ 


2.75(0.10) 


0.65 




0.65 


Wisp 1 


H 


3.15(0.20^ 
5.30(0.30|_ 


3.95(0.20) 


4.20 




6.00 


Wisp 2 


H 


6.25(0.25) 


8.40 




9.70 


Wisp 3 


H 


4.20(0.25) 


4.75(0.20) 


13.60 




14.00 


Anvil Wisp 1 


H 


3.20(0.15) 


3.60(0.15) 


5.30 




6.30 


Knot 


K s 


2.80(0.15) 


2.90(0.10) 


0.65 




0.65 


Wisp 1 


K s 


3.50(0.25) 


4.45(0.20) 


4.20 




6.00 


Wisp 2 


K s 


5.55(0.30) 


6.45(0.25) 


8.40 




9.70 


Wisp 3 


K s 


4.50(0.25) 


4.90(0.20) 


13.60 




14.00 


Anvil Wisp 1 


K s 


3.35(0.20) 


4.10(0.15) 


5.30 




6.30 
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Fig. 5. Spectra of wisp 2, the inner knot, together with spectra 
from the pulsar and the synchrotron nebula (Interwisp 1) for 
comparison. The pulsar spectrum (represented by X) is taken 
from Sollerman et al. (2002). The two infrared points in the spec- 
trum of the pulsar are from our infrared observations from 13 
December. The spectrum of the knot is shown by squares, the 
spectrum of wisp 2 is represented by diamonds, and the syn- 
chrotron nebula (Interwisp 1) is represented by triangles. All 
the data have been corrected for extinction using R v = 3.1 and 
E(B -V) = 0.52. 

ues of a v = -0.57 near the pulsar and decreasing to af v = -1 
at the edge of the synchrotron nebula. Veron-Cetty & Woltjer 
(1993) report values for a v for the region of the Crab PWN that 
we focussed on in the range -0.60 and -0.57, which is in good 
agreement with our results and which indicates that the spectrum 
has not changed since their measurements. 

4. Discussion 

We used continuum filters (U— >K S ) to check the emissivity of the 
Crab and then compared our results with the current theoretical 
models. The models from Del Zanna et al. (2003;2004;2006) 
seem to agree with our results. Del Zanna et al. (2006) man- 
age to recreate the observed features seen in the Crab nebula in 
both X-rays and the optical, and their modelled spectrum of the 
PWN shows a behaviour similar to our measured spectral index 
of the Crab PWN. Apart from the classical MHD approach, Del 
Zanna et al. (2006) added an additional equation that describes 
the evolution of the maximum energy of the emitting particles. 
This takes into account adiabatic and synchrotron losses along 
streamlines for the distribution of emitting particles. The spec- 



tral indices that they managed to estimate vary between -1.15 
and -0.6, the variation depending on the effective magnetisation 
of the wind plasma. 

We find that the structures within the Crab PWN (anvil and 
wisps) display spectra with similar spectral indices, although 
they are at the flatter end of the interval argued for by Del Zanna 
et al. Also, the regions between the wisps, despite being dim- 
mer, also display similar red spectra. The latter suggests that the 
wisps exist thank to local enhancements of synchrotron emis- 
sion. Two effects can be responsible for this. The first is that the 
wisps can be the result of an ion dominated wind that suffers 
a relativistic cyclotron instability (Spitkovsky & Arons 2004) at 
the termination shock. The other case is that the wisps are the re- 
sult of synchrotron cooling instabilities in a flow that is undergo- 
ing a transition from being particle dominated to field dominated 
(Hester et al. 1998;2002). 



4.1. Comparison with other pulsar wind nebulae 

There are PWNe around other pulsars that have been studied 
extensively. Data from Shibanov et al. (2003) for the Vela PWN, 
from Serafimovich et al. (2004) for the PWN of PSR B0540- 
69.3, and from Shibanov et al. (2008) for the PWN of 3C58 can 
be compared with the results for the Crab PWN. Table 6 shows 
the spectral indices we measured for the selected areas of the 
Crab PWN, along with measurements for the other PWNe. 

Shibanov et al. (2003) analysed images of the Vela pulsar and 
its surroundings in the J s and H bands and detected structures 
around the pulsar. One structure overlaps the inner arc region 
seen in X-rays (Helfand et al. 2001 ; Pavlov et al. 2001a) and can 
be the IR counterpart of the X-ray emitting region. No optical 
confirmation has been mentioned so far, but this can be due to the 
redness of this region (in the same manner as in the Crab), which 
makes it difficult to detect. The spectral index agrees with an 
extension from X-rays, i.e., a v in the range -0.5 to -0.3, which 
is marginally flatter than the spectra we derive for the inner part 
of the Crab PWN. 

A comparison between the Crab PWN and the PWN of PSR 
B0540-63.9 is of particular interest since this pulsar is often re- 
ferred to as the Crab twin due to its properties. Using HST im- 
ages from U to I, Serafimovich et al. (2004) measured the spec- 
tral index of the PWN around PSR B0540-63.9 to be a v « -1.48 
as a whole, while for various regions around the pulsar the spec- 
tral index of the PWN varies between a v = -0.27^ ^ and 
-1.5 8 ^ 43. The variation in the spectral index within the PWN 
of PSR B0540-63.9 depends on position (see Serafimovich et al. 
2004 for orientation), with a tendency for the brightest areas in 
the optical to have the steepest spectra. 
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Table 5. Photometric results from the Crab Wisps and the knot 
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Filter 


Inter Wisp 1 


Wisp 1 


Inter Wisp 2 


Wisp 2 


Inter Wisp 3 


Wisp 3 


Anvil Wisp 1 


Knot 


U 


1.70(0.10)^ 


1.70(0.10^ 


1.55(0.10) 


2.35(0.10)^ 


1.70(0.10) 


2.05(0.10) 


1.45(0.10) 




V 


1.70(0.10)" 


1.80(0.10) 


1.80(0.10) 


2.65(0.10) 


2.10(0.10) 


2.30(0.10) 


1.65(0.10) 




I 


1.75(0.10) 


1.95(0.10) 


1.85(0.10) 


2.85(0.10) 


2.25(0.10) 


2.40(0.10) 


1.75(0.15) 


1.60(0.10) 


H 


3.25(0.15) 


3.95(0.15) 


3.25(0.15) 


6.25(0.25) 


3.95(0.15) 


4.75(0.20) 


3.60(0.15) 


2.75(0.10) 


K s 


4.05(0.15) 


4.45(0.15) 


3.70(0.15) 


6.45(0.25) 


4.05(0.15) 


4.90(0.20) 


4.10(0.15) 


2.90(0.10) 



It should be kept in mind that, when comparing with the 
PWN of the distant PSR B0540-69.3 (for which 1" on the sky 
corresponds to the same true nebular size as -25" for the Crab), 
our Crab data probe regions much closer to the pulsar. For PSR 
B0540-69.3, a comparison should therefore also be made with 
the larger Crab survey of Veron-Cetty & Woltjer (1993), as dis- 
cussed and shown in Serafimovich et al. (2004). Nevertheless, 
it appears clear that the spectral index of the PWN around PSR 
B0540-69.3 is substantially steeper than that of the Crab PWN. 

3C58 a supernova remnant similar to the Crab (e.g., ple- 
rion) has been studied from radio (Weiler & Seiestad 1971) to 
IR (Shibanov et al. 2008). Using optical photometric data (B & 
V) and archive data from radio to infrared, Shibanov et al. (2008) 
constructed a spectrum for the PWN and the torus of 3C58. Their 
calculations of the spectral index for the PWN depend on the as- 
sumed extinction, but its value lies between a v = -0.5 (for the 
highest extinction) down to a v = -1.2. This brackets the X-ray 
spectral index of a v - -0.88 + 0.08, but there is no simple way to 
join the X-ray/optical/IR emission with a single power law. The 
same is true for the PSR B0540-63.9 PWN (Serafimovich et al. 
2004). 

In summary, there seems to be a slight increase in steepness 
of the optical/NIR spectra of the PWNe of Vela, Crab and 3C58, 
especially if one also includes the Crab IR results of Temim et 
al. (2006), who find a flattening off of the Crab PWN spectrum 
to between a v = -0.8 and a v = -0.3 around 4//m. The PWN of 
PSR B0540-63.9 is different from the other three in that it has 
a distinctly steeper spectrum, perhaps with the exception of the 
brightest region called "Area 6" by Serafimovich et al. (2004). 
As argued by Melatos et al. (2005), the cooling time for opti- 
cal/NIR photons is long compared to the flow time for the inner 
part of the Crab PWN. Based on the similar spectral slopes, it 
is likely that the same could also be true for Vela and 3C58, 
whereas cooling could be more efficient for most parts of the 
PSR B0540-63.9 PWN (see also Serafimovich et al. 2004) than 
for the other PWNe. 



5. Conclusions 

We have presented data for the Crab pulsar in both the optical 
and infrared. Our infrared data suggest changes in the inner en- 
viroment of the Crab for both dynamic structures and emissivity. 
Our estimated velocities for the wisps do not agree with cur- 
rently known values (Hester et al. 2002). The reason for this 
difference can be the large time window we had. Ideally, the 
dynamic structure of the inner Crab should be studied using a 
daily based program. From our data we find that the inner knot 
does not show any variability in its emission and that it remains 
stationary over a period of 13.5 years (until 2007 December) rel- 
ative to the pulsar; its average projected velocity relative to the 
pulsar within this time window is < 8 km s" 1 , corresponding to 
a relative space motion of < 10 km s" 1 , if the pulsar and knot 
are aligned with the jet axis. 



Table 6. Measued spectral indices from known PWN 



Pulsar (Area) 


Wavelength range 


Spectral Index 


Crab (Wisp 1) 


K s 


-0.58 + 0.08^ 


Crab (Wisp 2) 


K s 


-0.63 + 0.07 


Crab (Wisp 3) 


K s 


-0.53 + 0.08 


Crab (Anvil Wisp) 


K s 


-0.62 + 0.10 


Crab (Interwisp 1) 


K s 


-0.61 + 0.08 


Crab (Interwisp 2) 


K s 


-0.50 + 0.10 


Crab (Interwisp 3) 


K s 


-0.49 + 0.10 


Crab (Knot) 


I-> K s 


-0.63 + 0.02 


Vela (Inner arc)^ 


X-rays^IR 


-0.5 -> -0.3 


3C58 (Optical region)^ 


X-rays 


-0.88 + 0.08 


3C58 


Optical^IR 


-1.2 ->-0.5^ 


PSR B0540-69.3 (Area 1)^ 


U->I 


_i 09+ - 30 


PSR B0540-69.3 (Area 2) 


U->I 


i CO +0.42 


PSR B0540-69.3 (Area 3) 


U->I 


-1.28*$ 


PSR B0540-69.3 (Area 4) 


U->I 


_i 4Q+0.59 


PSR B0540-69.3 (Area 5) 


U->I 


-0 Q2 +0 - 33 


PSR B0540-69.3 (Area 6) 


U->I 


-0 27 +0 9 ° 

U ' Z/ -0.91 



The combination of the optical and NIR data allowed us to 
create spectra of the wisps and the knot and thus obtain hints to 
their nature. All the wisps (and the interwisp medium) display 
red spectra with similar indices (a v that vary between -0.49 ± 
0.10 and -0.63 ± 0.07), and the knot has a v = -0.63 ± 0.02. 
The latter seems to agree with the result of Sollerman (2003), 
who reported a spectral index for the knot equal to a v = -0.8 
for a spectral range that covers both optical (U, V and I) from 
the HST and IR (J, H, K s ) from the VLT. The spectral distribu- 
tions we find for the wisps and the red knot are similar to what 
has been estimated for other PWNe, which suggests that these 
observed properties are the result of the same emission mecha- 
nism. The exception appears to be the PWN of PSR B0540-63.9, 
which for most of its PWN has a steeper spectrum than the other 
PWNe and which signals stronger cooling. A feature in the Crab 
PWN that sticks out is the red knot with its high degree of po- 
larisation. It lies close to the pulsar and could be part of a shock 
structure shared with the sprite. Both the pulsar and the sprite 
have markedly flatter spectra than the knot though. On a some- 
what larger scale, we find that our data are compatible with the 
model of Del Zanna et al. (2006), although our spectra are at the 
flatter end of the spectral index interval argued for by Del Zanna 
et al. The similarity of the spectral indices of the wisps with the 
interwisp medium can be an indication that the wisps are just 
enhancements of synchrotron emission, but synchrotron cooling 
instabilities should not be dismissed. 

We highlighted that pulsar wind nebulae of Crab like rem- 
nants display similar red spectra, and observations have also re- 
vealed that in some cases they show similar structures to those 
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of the Crab (e.g., Vela in X-rays). Further detailed comparisons 
are needed to constrain the physical processes giving rise to the 
observed properties. Because there are many open questions re- 
garding the nature of the PWNe, it would be interesting to extend 
this comparison to related objects like Crab-like remnants with 
bow- shock pulsars such as the Guitar Nebula, in particular to 
study how age affects the properties of the pulsar wind nebulae. 
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